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Abstract 

We notice that, when a quantum system involves exceptional points, i.e. the special 
values of parameters where the Hamiltonian loses its self-adjointness and acquires 
the Jordan block structure, the corresponding classical system also exhibits a sin- 
gular behaviour associated with restructuring of classical trajectories. The system 
with the crypto-Hermitian Hamiltonian H = (p 2 + z 2 )/2 — igz 5 and hyper-elliptic 
classical dynamics is studied in details. Analogies with supersymmetric Yang-Mills 
dynamics are elucidated. 

1 Introduction 

A sufficient condition for the spectrum of a Hamiltonian to be real is its self-adjointness, 
= H . However, there exists a rich class of systems whose Hamiltonian is not manifestly 
Hermitian, but the spectrum is still real. Probably, the simplest such example is the 
matrix Hamiltonian 

with two real eigenvalues Ai = 1 and A2 = 2. Such systems have been intermittenly 
discussed since mid-seventies [1] , but the interest to this problem was considerably boosted 
after the paper [2] appeared, where the systems with complex potentials V(x) = x 2 (ix) e 
were studied and it was shown that the spectrum of the corresponding Hamiltonians is 
real. 

Generically, a Hamiltonian involving only real eigenvalues can be transformed in a 
manifestly Hermitian form by a similarity transformation [3], H — > RHR^ 1 . If H is 
not manifestly Hermitian, R is not unitary. This amounts to modifying the Hilbert 
space measure so that the probability P = [ip\M\ip) defined with the new measure M is 
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conserved and the theory is unitary. One can call a Hamiltonian of this type crypto- 
Hermitian (Hermitian in disguise) [4,5]. 

However, not all Hamiltonians with real eigenvalues are crypto-Hermitian. In some 
cases, a Hamiltonian cannot be rendered Hermitian (and eventually diagonalized) by a 
similarity transformation. The simplest example is a 2 x 2 matrix representing a Jordan 

block, H — ( ^ ). The exceptional points when this happens are associated with degen- 
eracy of eigenvalues [6]. Exceptional points have measure zero in the space of parameters. 

There are exceptional points involving multidimensional Jordan blocks. For the higher- 
derivative Pais-Uhlenbeck oscillator [7] at equal frequencies, the dimension of emerging 
Jordan blocks is even infinite [8, 9] . A more typical kind of exceptional point is when 
only a couple of eigenvalues coalesce. Changing the parameters beyond the exceptional 
point, a pair of complex eigenvalues that are conjugate to each other appears. @ 

A numerical study performed in Ref. [2] displayed an infinite number of exceptional 
points in the parameter e for the Hamiltonian 

H = p 2 + z 2 (tzY . (3) 

The exceptional points lie in the interval e G (—1,0). The leftmost exceptional point 
—.578. When e < e*, the spectrum of the Hamiltonian involves only one real 
eigenvalue (the ground state), with all other eigenvalues coming in complex conjugate 
pairs. At e = e*, one of such pairs coalecses and, at still larger e, is transformed into a 
pair of close real eigenvalues. Then comes the turn of the second pair, etc. The infinite 
number of exceptional points has an accumulation point at e = [10]. At e > 0, complex 
eigenvalues disappear and the Hamiltonian is crypto-Hermitian. 

A similar phenomenon was observed in Ref. [11] for the Hamiltonian 



There are two different spectral problems corresponding to this Hamiltonian (see Refs. 
[11, 12] for detailed explanations). One of the problems is defined in the Stokes wedges 
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1 The explicit forms of R and M for the Hamiltonian |l| are 

*= (J I') > M = BFR= (_\ - 1 ) (2) 

2 Again, the simplest example is the matrix H = ( | ). When a is small and positive, we have 

a pair of close real eigenvalues. When a is small and negative, there is a pair of complex conjugate 
eigenvalues, a = is the exceptional point. 
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including the real axis. It turns out that the spectrum is real there for all values of g. 
But for another spectral problem defined in the wedges 

57T , . 7T 

- n < Arg(z) <- n , 

157T . . 197T . 

— < Arg(z) <— (6) 

in the complex z plane, the situation is different. The spectrum is real (and the Hamilto- 
nian is crypto-Hermitian) for large enough g, but the system involves an infinite number 
of exceptional points at small g. If g < g* ~ 0.037, a pair of complex eigenvalues appears. 
If g < g„ ~ 0.007, there are at least two such pairs, etc. 

A question that can be asked is whether the presence of these quantum exceptional 
points displays themselves in some way also in the dynamics of the corresponding classical 
systems ? A partial answer to this question was obtained in Ref. [13] where a series of 
critical values of the parameter e, where the pattern of classical trajectories in the problem 
(J3]) is changed, was found. The values of the classical exceptional points do not coincide 
with the values of the quantum exceptional points. One can only say that the presence 
of the former is associated with the presence of the latter. 

The problem is rather complicated, however. The Riemann surface for the potential 
~ z 2+€ has generically an infinite number of sheets, the classical trajectories can visit a 
lot of the sheets and look complicated. That is why we decided to study this question for 
the system (j3J), which is much simpler. The classical trajectories represent in this case 
hyper-elliptic functions known to mathematicians [14]. 

The result is the following. There are no classical exceptional points for the system (J3J) 
with a positive energy. When the energy is negative, there is a single classical exceptional 
point 
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As was mentioned above, a classical exceptional point is the point where the pattern of 
the classical trajectories changes. For the system (j3J), the reason for this change is very 
simple: it happens that at g — gf^ ss two of the five turning points, i.e. two of the five 
solutions to the equation 

z 2 

V(x) = igz 5 = E < (8) 

coalesce. This phenomenon has a lot in common with the so called Argyres- Douglas 
phenomenon in supersymmetric Yang-Mills theory [15]). We will dwell upon this issue in 
Sect. 4. 

Before discussing the dynamics of the system (T4j), we want to make few comments on 
the classical dynamics of a simpler system 

H-tlgt+tf, (9) 

where exceptional points do not appear either at the classical or at the quantum level. 
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2 Complex cubic potential and elliptic dynamics. 



Consider the system Complex classical trajectories are the solutions to the equation 

z 2 UJ 2 Z 2 o . s 

y + ~y~ + ^ 3 = • ( 10 ) 

The solutions to this equation with real positive energies (especially, in the case u = 0) 
were found numerically in Ref. [2]. The solutions with real negative energies were found 
in Ref. [11]. A simple remark that we want to make here is that the solution to Eq. ffTOl) 
has a name - it is the Weierstrass function. Let for simplicity set u — 0, g — 1. The 
equation is brought into the canonical Weierstrass form [16], 

y 2 = Ay 3 - g 2 y - g 3 , (11) 

with g 2 = 0, g 3 = E/2, by setting z = 2iy. 

The Weierstrass function z = P{t,g 2 ,g 3 ) is an elliptic double-periodic complex func- 
tion. Generically, the periods are complex, but in our case one of the periods is real. For 
positive energies, 

Ti = 5 4 mjk E ~ lle K 3AME ~ 1 " ■ (12) 

Another period is still complex, T 2 = Tie" 1 "/ 3 . T 2 is obtained from 7\ by rotation in 
the complex E plane, E — > Ee~ 2t7T . Any linear combination of T\ and T 2 with integer 
coefficients is also a period. A purely imaginary combination iT = 2T 2 — T\ — %\[ZT\ 
represents a particular interest. The physical interpretation of the periods T\ and T is 
clear [2, 11]. The periodicity with respect to the real time shift t — > t + T\ is a usual 
physical periodicity of the trajectories. For all trajectories, the period is the same. The 
imaginary time shift transforms one trajectory into another. [f| 

One can also be interested with the solutions to Eq. (1101) at complex E. It is still 
Weierstrass function, but all its periods are now complex. This means that the solution 
is not periodic in real physical time. Such solutions were studied numerically in Ref. [17]. 

The family of all trajectories contains two distinguished members that can be called 
stem trajectories. One of the stem trajectories connects two turning points (solutions 
to the equation V(z) = iz 3 = E) in the lower half-plane, Z\ = E 1 / 3 e~ 5l7T / 6 and z 2 = 
jr;i/3 e -in/6 _ Another stem trajectory goes from the turning point z 3 = iE 1 ^ 3 to infinity. 
The shift t — > t+iT/4 transforms one stem trajectory into another. The shift t — > t+iT/2 
transforms the function z = 2iP(t, 0, E/2) to z = 21P(— t, 0, E/2) and the shift by iT 
leaves the Weierstrass function intact. This is all illustrated in Figs. 1,2. 

Translating these observations into a standard mathematical language, one can say 
that the Weierstrass function describes the motion over the Riemann surface of the func- 
tion \[E — iz 3 . The turning points and infinity are nothing but the branching points 

3 This transformation can be interpreted as a gauge transformation [11], but it is not a point of interest 
for us now. 
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Figure 1: Periods of the Weierstrass function. 



Figure 2: The trajectories z = 2iP(t + ia, 0, 1/2), t £ (0, T\). The vertical line (a = 0) is 
the stem trajectory connecting the branching points z = i and z = oo. The "smile" (a = 
T/4) is the stem trajectory connecting the branching points z = e~ 5i7r / 6 and z = e~ l7T / 6 . 
The trajectories with intermediate values a = T/16 and a = T/8 are also plotted. 

of this function. This Riemann surface has the topology of torus. Two periods T 12 
corresponds to two cycles of this torus. 

One can also plot the imaginary time trajectories 2iP(a + it,0, 1/2), t £ (0,T) and 
observe by inspecting Eq. lflTT) that they coincide up to a sign with the real time trajectories 
2iP(t — ia, 0, —1/2) having negative energies and studied in Ref. [11]. 

3 The potential V(z) = z 2 /2 — igz 5 and hyper-elliptic 
dynamics. 

Classical complex trajectories for the Hamiltonian have much in common with the 
elliptic trajectories of the previous section. They are hyper-elliptic trajectories decribing 
the motion over the Riemann surface of the function E — V(z). This Riemann surface 
has genus 2 and two pairs of cycles. The latter implies that the equation 

Y + V(z) = E (13) 
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has now different types of solutions stemming from the trajectories that connect three 
different pairs of branching points of the function yf E — Vjzj. Let first the energy be 



positive and V(z) 



-iz 



(without the quadratic term). The branching points are in the 



vertices of the pentagon and at infinity. Three different stem trajectories are depicted in 
Fig.3. The corresponding periods are [2, 11] 



Figure 3: Stem trajectories for the potential —iz 5 
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and T 3 = T\ — T 2 . 

Similarly to what we had for the cubic potential, the solutions to Eq. (fT3|) are also 
periodic with respect to imaginary time shifts, 
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The periods (114)15p are all interrelated under rotations in the complex E plane (the 
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monodromies). Q For example, 



Tl (Ee~ 2 n = T 2 (Ee- 2 n = (16) 

This leads to three families of trajectories with positive energies and three families of 
trajectories with negative energies. U 

Let us switch on the quadratic term now and investigate how the pattern of trajectories 
is changed when changing g. Let first the energy be positive. Note that for all g, Eq.(jHJ) 
still has 5 distinct roots (see Fig. 3). 
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Figure 4: Turning points. E = 1, g = 0.1 



We see that the turning points are still symmetric with respect to the imaginary axis. 
There are still three families of trajectories stemming from the trajectories that have the 
same qualitative pattern as in Fig. [3j Nothing changes essentially when g is changed. 

When E < 0, the situation is different. There is a point (j7|) where two of the roots 
of Eq.(jSJ) coalecse. The root patterns slightly above the exceptional point g = g* and 
slightly below it is shown in Fig. [5j 

When g > g*, there are three families of orbits. Their stem trajectories are shown in 
Figj6h.. When g < g*, the pattern of the stem trajectories is different (see Fig. [6b). One 
can notice that the stem trajectory connecting the low pair of the turning points with 
Re(x) 7^ 0, exists at both g > g* and g < g*, but its form is different. When g — > g* + iO, 
the trajectory crosses itself, after which the "appendix" going around the low turning 
point at the imaginary axis disappears. 

4 Cf. Proposition 7 in Ref. [18]. 

5 To be more precise, the third family is degenerate in this case and consists in only one member: the 
stem trajectory starting at z = iE 1 ' 5 and going away to z = ioo in finite time, < C scapc = T3/2. 
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a) b) 
Figure 5: Turning points. E = — 1; a) g = 0.06, b) g = 0.03. 

4 Analogies with SYM. 

As was already mentioned, the restructuring of classical trajectories observed in this paper 
has much in common with the Argyres-Douglas phenomenon that occurs in SYM theory. 
Let us briefly explain this point here (see Ref. [19] for more details). 

M = 2 supersymmetric gauge theories are known to include vacuum moduli spaces: 
valleys of degenerate classical vacua parametrized by vacuum expectation values of scalar 
fields [20]. These valleys have some distinguished singular points where the mass of 
certain objects present in the theory (the monopoles and dyons) vanishes. The original 
AD observation was that, in some theories and for certain values of parameters, the moduli 
space singularities might coalesce. 

An analogy between this phenomenon and the classical exceptional point ([7]) exists, 
but is not so manifest. It becomes much more clear for the theories where J\f = 2 is 
broken down to Af = 1. The continuous degeneracy of the vacuum valleys is lifted in 
this case, and we are left with a finite number of different classical vacua. These vacua 
can be separated by domain walls, the classical solutions to the equations of motion with 
proper boundary conditions. Mathematically, domain walls interpolating between the 
classical vacua play exactly the same role as the stem trajectories interpolating between 
the turning points in the QM models discussed above. 

The point is that, in some M = 1 models, the parameters can be chosen such that 
the classical vacua coalesce. This phenomenon is akin to the original AD phenomenon 
(merging of singularities in vacuum moduli spaces). One can call it the AD phenomenon 
of the second kind. It was discovered first for a somewhat exotic model based on the 
G2 gauge group and involving besides J\f = 1 gauge supermultiplet three different chiral 
multiplets S^"l 1 1 ' 2 ' 3 7 in the fundamental representation of G2 [21]. The superpotential of 
the model is 

777 A 

W = -y(^) 2 - ^f^S^S* , (17) 
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a) b) 

Figure 6: Stem trajectories. E = — 1; a) g = 0.06, b) g — 0.01. The small gaps between 
the vertical lines going to ioo and the curves below are indistinguishable. 

where f" 13 " 1 is the invariant antisymmetric tensor of G2 u- To analyze vacuum struc- 
ture, one has to add to the superpotential (TIT]) a nonperturbative instanton-generated 
superpotential of the Affleck-Dine-Seiberg type [22]. In this case, it has the form 

A 9 

£ 2 -det.M' 1 ' 

where B and M. are the gauge-invariant moduli, 



M ij = S* a Si, B = \^ k S l a S^ (19) 



and A is the confinement scale. 

Adding the superpotentials ffTTI) and (fTSj) . we obtain a theory with two dimensionless 
parameters: the Yukawa coupling A and the ratio mass m/A. The equation determining 
the chirally asymmetric classical vacua is of the 6-th order, 

SW 0^W(1-^V = A°, (20) 



^(moduli) \ m? 

where u = |TrA^. Generically, there are 6 distinct solutions. However, when A 2 ~ 
.385m 2 /A 2 , two of these vacua coalesce. 

The simplest model where the AD phenomenon of the second kind occurs is based on 
SU(2) gauge group and involves a massive adjoint chiral multiplet and a pair of funda- 
mental chiral multiplets. This model obtained by a deformation of Af = 2 supersymmetric 



6 Note that G2 can be defined as a subgroup of 50(7) leaving invariant the form f a ^^ ' A a BfiC 1 for any 
set of three different 7- vectors A,B,C. 

7 We leave aside a controversial issue of possible existence of the extra chirally symmetric vacuum [23] . 
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QCD (the model involving besides the gauge Af = 2 multiplet also a matter hypermulti- 
plet) and analyzed in Ref. [24] involves generically 3 vacua, but, for some values of masses 
and Yukawa couplings, two of them coalesce. 

When the vacua merge, the domain wall interpolating between them disappears. The 
analogy with the classical exceptional point d7j) is obvious. 

I aknowledge numerous illuminating discussions with A. Gorsky and warm hospitality 
at AEI institute in Golm where this work was finished. 
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